Collective modes and ballistic expansion of a Fermi gas in the BCS-BEC crossover 
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We evaluate the frequencies of collective modes and the anisotropic expansion rate of a harmon- 
ically trapped Fermi superfluid at varying coupling strengths across a Feshbach resonance driving 
a BCS-BEC crossover. The equations of motion for the superfluid are obtained from a microscopic 
mean-field expression for the compressibility and are solved within a scaling ansatz. Our results 
confirm non-monotonic behavior in the crossover region and are in quantitative agreement with 
current measurements of the transverse breathing mode by Kinast et al. [Phys. Rev. Lett. 92, 
150402 (2004)] and of the axial breathing mode by Bartenstein et al. [Phys. Rev. Lett. 92, 203201 
(2004)]. 
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Current experiments on ultracold Fermi gases are 
rapidly advancing towards the realization of superfluid 
states, and Bose-Einstein condensation of dimers has al- 
ready been achieved P, 0> 0, 0. IE 01- A key tool for 
the manipulation of atomic gases is the use of a Fesh- 
bach resonance to vary the magnitude and sign of the 
coupling strength. Across the resonance the s-wave scat- 
tering length a goes from large positive to large negative 
values, thus allowing exploration of the crossover from 
the Bardeen-Cooper-Schrieffer (BCS) state to the Bose- 
Einstein condensate (BEC) of bound-fermion pairs. As 
Fermi gases have been demonstrated to be stable also 
near the resonance 0, they offer a new opportunity to 
investigate highly correlated many-body systems. 

As in the case of bosonic clouds, the frequencies of col- 
lective modes of Fermi gases can be measured to high 
accuracy and can yield information on the state of the 
system and on its collisional properties, thus contribut- 
ing to its characterization in a strong-coupling regime. 
Another dynamical observable which is readily accessed 
in the experiments is the aspect ratio of an expanding 
cloud released from an anisotropic trap, which depends 
on its quantum state and on its degree of collisionality. 
While a normal Fermi gas in the collisionless regime ex- 
pands spherically, an anisotropic expansion is predicted 
for a normal gas in the collisional regime as well as for a 
superfluid gas. 

The purpose of this Letter is to present a theory of the 
collective modes and of the expansion of a trapped su- 
perfluid Fermi gas at zero temperature as the coupling 
strength is varied across a Feshbach resonance. The 
frequencies of the collective modes are well known in 
both the BCS and the BEC limit HUE!. From non- 
mean-field perturbative estimates it has been conjectured 
that the frequency of the transverse breathing mode in a 
highly elongated trap should exhibit a non-trivial depen- 
dence on the scattering length Two further studies 
have been based on semi-empirical forms of the equation 
of state 0, ^| . Here we use a microscopic mean-field 
description of the BCS-BEC crossover 0, which 



physics in all regimes ^(|. We calculate the equation 
of state and the density profiles of the gas under axially 
symmetric confinement with the help of a local density 
approximation (LDA), and use them to determine the 
collective mode frequencies and the expansion rate by 
means of a simple scaling assumption. Our results do 
not use an interpolation scheme nor involve adjustable 
parameters, and show already at mean-field level non- 
monotonic behaviors across a Feshbach resonance. 

The frequencies of radial and axial breathing modes 
have very recently been measured in elongated clouds of 
6 Li atoms at various values of the coupl ing strength in 
the strong-coupling intermediate regime 

[HEl. Quite 

remarkably, our mean-field approach yields predictions 
for the frequencies of the breathing modes that are in 
good quantitative agreement with these experiments. 
Equation of state and equilibrium density profile. 
Our starting point is to determine the chemical poten- 
tial fx (n) as a function of density n for a homogeneous 
Fermi gas through the BCS-BEC crossover at zero tem- 
perature. We use a mean-field theory proposed in the 
pioneering work of Leggett 0] and of Nozieres and 
Schmitt-Rink [l|. This describes well both the BCS 
weak-coupling regime and the BEC strong-coupling limit, 
and is believed to be a reliable approximation in between 
[Til Il5l ITtH . It has also been applied to a trapped gas 
within the LDA [l9ll20|. 

The gas consists of two equally populated spin com- 
ponents that interact by a contact pseudopotential 
parametrized by the s-wave scattering length a. The 
mean-field theory for the crossover extends the usual BCS 
gap equation 
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FIG. 1: (color online) The exponent 7 as a function of the 
dimensionless coupling parameter k = (N 1 ^ 6 a/aho)^ 1 for a 
trapped Fermi gas with N — 2 x 10 5 and A = lo z /lu± — 
0.05. The quantity 7 (n) and the transverse density profile 
no = no (p, 0) with p — x 2 + y' 2 are plotted in the insets for 
k = -1.0, +0.1, +0.5 and +1.0 (from bottom to top). The 
coordinate p is in units of auo — X^ 1 ^ 6 yj / h/mcu± and no is in 
units of a7^. 



tion this exponent is independent of the index a and it is 
seen from Fig. 1 that it decreases in going from the BEC 
to the BCS limit — i.e., the trapped gas is relatively less 
compressible in the BEC limit. 

Scaling ansatz. — We describe the dynamics of density 
fluctuations in a superfiuid Fermi gas by the equation of 
continuity 



d t n + V • (nv) = 0, 



(3) 



and the Euler equation 

mdtv + V [n (n) + V ext + mv 2 /2] = 0. (4) 

Here n (r, t) and v (r, t) are the time-dependent den- 
sity profile and velocity field of the trapped gas, [i (n) 
is the local chemical potential expressed as a func- 
tion of n (r, t) in an adiabatic LDA, and V ex t (r, t) — 
m [w 2 (t) (x 2 + y 2 ) + lu 2 (t) z 2 ] /2 is the axially symmet- 
ric trapping potential. The above equations apply as well 
to a normal Fermi gas in the collisional regime. The dif- 
ference from a superfiuid should emerge from some fur- 
ther knowledge of the collisional behavior of the cloud or 
from its response under rotation. 

We adopt the following scaling form of the time- 
dependent density profile, 



h 2 k 2 /2m, .Ek = \J (tk — A*) 2 + A 2 , A is the gap parame- 
ter, and n is the total density for fermions in either spin 
state (n-t- = n\ = n/2). These equations are to be solved 

for fi and A with a given choice of the dimensionless 

1/3 

coupling parameter kpa = (37r 2 n) a. In the weak- 
coupling limit (kpa — > CP) they give back the standard 
BCS result, and when fc^a — > + they correctly repro- 
duce the binding energy of fermion pairs at leading order 
[l6[ . In the unitarity limit (a = ±00) the model gives 
/j,(n) (x n 2 / 3 , scaling as predicted by the universality hy- 
pothesis |2ll |. 

We then use the equation of state fi(n) to deter- 
mine the density profile Hrj(r) under the confining po- 
tential V ex t(r) — m [lo^x 2 + y 2 ) + muj 2 z 2 ] /2 through 
the implicit LDA equation /i(no(r)) + V ex t(^) = fa- 
Here fa is the chemical potential of the trapped gas 
from the normalization condition N = J ' d 3 r tiq (r), N 
being the total number of fermionic atoms. A useful 
dimensionless coupling parameter for a trapped gas is 
k = (N^a/aho)- 1 pb 1.695(fc» 
Fermi wave number at the trap center [TT| 

Our results are illustrated in Fig. 1. The left inset re- 
ports the function 7(77.) = 1+n [(cPfi/dn 2 ^) / (dfi/dn)] in 
the homogeneous gas for four choices of k, while the right 
inset shows the LDA density profiles in the transverse di- 
rection. In the special case of a power-law dependence of 
the chemical potential on density n(n) oc rt 7 , j(n) = 7. 
In the general situation 7(71) acts as an effective expo- 
nent. The main body of Fig. 1 shows the quantity 
7 = (/ d 3 rn r 2 7 ("0)) / (/ d 3 rn r 2 ) with a = x, y, or 
z, to be used below as an averaged effective exponent for 
inhomogeneous gas. In fact, in the present approxima- 
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where hq (r) is again the equilibrium density profile and 
the dependence of the profile n (r, t) on time is entirely 
contained in the scaling parameters b a (t) [22|. This 
ansatz correctly reflects the expansion and compression 
along each axis and therefore is appropriate for the de- 
scription of the breathing modes and of the ballistic ex- 
pansion of the cloud. The corresponding form of the 
velocity field is fixed by the equation of continuity as 
v a (r,t) = b a (t)r a /b a (t). 

The scaling ansatz in Eq. J5J is an exact solution of 
the equations of motion if the equation of state is a power 
law pi (n) oc n 7 . In this case the scaling parameters obey 
the coupled differential equations 



b a + cj 2 b a - (u>l/b a 




= 0. 



(6) 



However, with a more general form of the equation of 
state the scaling solution is not satisfied at every posi- 
tion r. A useful approximation is to take a spatial av- 
erage with the weight function r 2 no [2i|. After some 
straightforward algebra we obtain 
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FIG. 2: (color online) The frequency of the axial and trans- 
verse breathing modes in an elongated trapped Fermi gas as 
functions of k = (N 1/ ' 6 a/af lo )~ 1 . The solid circles and empty 
triangles with error bars are the experimental results given by 
Kinast et al. [l7j | and by Bartenstein et al. [Tsl ]. respectively. 
Here we plot only the measured frequencies with low dacay 
rates in the resonance region, for instance, for the axial mode 
we use the data of the inset of figure 1 in Ref. flSj, The 
resonance position Bo has been set at 837 G. The error bar 
in the experimental data in the horizontal (re) axis is due to 
the uncertainty of Bo and we have taken A£>o = ± 15G. 



which reduces to Eq. (j^J if n (n) oc n 1 . 
Breathing modes — We evaluate the transverse and ax- 
ial breathing modes by taking periodic variations of the 
trapping frequencies u>± ;Z (t) = lu±. z (1 + ecos(wi)) with 
e <C 1. Linearization of Eq. (0 yields the linearized 
form of Eq. © with 7 replaced by 7. One immedi- 
ately finds three mode frequencies, one of which lies at 
lj = \/2uj^ independently of the coupling strength and 
belongs to a surface mode with projected angular mo- 
mentum m — ±2. The other two mode frequencies are 
given by 



: >1 = l+7 + (2 + 7 )A 2 /2 



± V [1 + 7 + (2 + 7)A 2 /2] Z - 2 (2 + 3 7 ) 

where A = uj z /uj± and the ± signs refer to the transverse 
and axial mode, respectively. 

In the highly elongated traps of current experimen- 
tal interest ( A <C 1) the frequencies of the breath- 
ing modes reduce to u> + /tu± = y/2 + 27 and ui-/u> z = 
yj{2 + 37)/ (I + 7), yielding uo + = ^/3uj ± ~ 1.83w_l 
in a dilute Fermi gas (7 = 2/3) and lo + = 2uj±_ in a dilute 
BEC (7 = 1). The mode frequencies in the crossover 
region are presented in Fig. 2 as solid lines displaying 
a non-trivial dependence on the coupling parameter. In 
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FIG. 3: (color online) The frequency of the transverse and ax- 
ial breathing modes as functions of the anisotropy parameter 
A for a trapped Fermi gas with N = 2 x 10 5 , at the indicated 
values of the coupling constant re. 



particular, a dip appears near the unitarity limit on the 
BCS side, as a result of the pairing that enhances the 
compressibility of the gas. On the BEC side the frequen- 
cies monotonically increase towards the BEC value, this 
behavior being related in our approximate scheme to the 
increase of the effective exponent 7 with the coupling con- 
stant k for k > 0. In Fig. 2 we have also compared our 
predictions for the axial mode with experimental data 
by Bartenstein et al. 0] and those for the transverse 
mode with measurements by Kinast et al. |l7j| and by 
Bartenstein et al. [lq. Two sets of the experimental 
data, i.e., the axial mode in Ref. and the transverse 
mode in Ref. [l7| , are in quantitative agreement with our 
mean-field results. 

However, the remaining set of data for the trans- 
verse mode measured by Bartenstein et al. con- 
tradicts our prediction. The measured frequencies have 
an abrupt change around the resonance, accompanied by 
large damping rates. This is not expected by the the- 
ory. As noted by Grimm, this might be due to the ac- 
cidental coincidence between the collective energy and 
the low-lying quasiparticle energy. As a result, the su- 
perfluidity of the Fermi gas could be destroyed by the 
excitations of quasiparticles. Nevertheless, the tendency 
of monotonic increase of the predicted mode frequency 
in the BEC limit, as a natural outcome of the mean-field 
theory, agrees qualitatively with the experimental obser- 
vations. This is to be contrasted with the findings by 
Stringari which predict a decrease in the transverse 
mode frequency towards the BEC limit. 

In our calculations the collective mode frequencies in 
the trap depend mainly on the dependence of the chem- 
ical potential on density through an effective power-law 
exponent. We have not taken account of non-mean-field 
corrections to the compressibility, which were discussed 
in Ref. ^lj. These can be estimated in a perturbative 
way for re — > ±00, and can be shown to be very small 
near the BEC limit from the Lee-Huang- Yang correction 
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FIG. 4: (color online) The aspect ratio Xb±/b z as a function 
of k at a fixed time r = ui±t — 25. The inset shows the aspect 
ratio as a function of the dimensionless time r for four values 
of the coupling parameters k as indicated. For comparison 
the expansion of an ideal gas is also shown in the inset by a 
thin solid line. The other parameters are N = 2 X 10 and 
A = 0.05. 

term [24J evaluated with an estimate of the molecule- 
molecule coupling constant by Petrov et al. 25]. The 
Hartree term, on the other hand, tends to emphasize the 
dip shown in Fig. 2 near the BCS limit. 

The calculated collective modes as functions of the 
anisotropy parameter A are shown in Fig. 3 for four val- 
ues of the coupling strength. The sensitivity of the trans- 
verse mode to the coupling is preserved for all values of A, 
while the axial mode is less sensitive to the interactions. 
For a spherical trap with (A = 1) the breathing modes are 
known as the monopole and the quadrupole mode and we 
have LO+/u± = \/1 + 3j and u^/uj z = y/2. In the oblate 



limit (A > 1) the frequencies are w+/w z = -y/2 + 7 and 
u>_/u ± = v /(4 + 6 7 )/(2 + 7 ). 

Expansion of the Fermi gas. - Finally we turn to 
study the expansion of the Fermi gas after suddenly 
switching off the confining potential, by setting u) a (t > 
0) = 0. We follow the time evolution of the aspect ra- 
tio R±(t) I R z {t) = Xb± (t) /b z (t) by numerically solving 
Eq. (J2J with the initial configurations b a (0) = 1 and 
&«(0) - 0. 

The results are shown in Fig. 4. The aspect ratio of the 
cloud at a given expansion time shows non-monotonic be- 
havior as a function of the coupling strength k, similarly 
to what is observed in the collective modes frequencies. 
This reflects again the emergence of important interac- 
tion effects in the hydrodynamic behavior. 
Conclusions. — In summary, we have presented a micro- 
scopic theory of the dynamics of a superfluid Fermi gas 
through the whole BEC-BCS crossover. Our analysis is 
based on the equation of state of the gas in a mean- 
field approach. For a gas under anisotropic confinement 
we have found that both the frequencies of the breath- 
ing modes and the aspect ratio of the expanding cloud 
show non-monotonic behavior as functions of the cou- 
pling constant. We have also made quantitative contact 
with current experiments. 
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